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Abstract. We prove that the number of Hamilton cycles in the random graph 
G(n,p) is n!p n (l + o(l)) n a.a.s., provided that p > lnn + lnl ^ 1 "+"( 1 ) _ Furthermore, we 
prove the hitting-time version of this statement, showing that in the random graph 
process, the edge that creates a graph of minimum degree 2 creates (^p)™ (l + o(l))" 
Hamilton cycles a.a.s. 



1. Introduction 

The goal of this paper is to estimate the number of Hamilton cycles in the random 
graph G(n,p). To be more formal, we show that the number of Hamilton cycle is a.a.s. 
concentrated around the expectation up to a factor (1 + o(l)) n , provided the minimum 
degree is at least 2. 

It is well known (see e.g. [5]) that the minimum degree of G(n,p) is a.a.s. at most 
one for p < lnn + lnln "~^( 1 ) ; anc [ therefore G(n,p) contains no Hamilton cycle in this 
range of p a.a.s. Komlos and Szemeredi [23] and Korshunov [24] were the first to show 
that this bound is tight, i.e., G(n,p) is a.a.s. Hamiltonian for every p > lnn + lnlnn + a; ( 1 ) _ 
Bollobas [6] and independently Ajtai, Komlos, and Szemeredi [1] proved the hitting 
time version of the above statement, showing that in the random graph process, the 
very edge that increases the minimum degree to two also makes the graph Hamiltonian 
a.a.s. 

There exists a rich literature about hamiltonicity of G(n,p) in the range when it is 
a.a.s. Hamiltonian. Recent results include packing and covering problems (see e.g. [13], 
[21], [22], [26], [15], and [18]), local resilience (see e.g. [34], [14], [4] and [28]) and Maker- 
Breaker games ([33], [16], [3], and [11]). In this paper, we are interested in estimating 
the typical number of Hamilton cycles in a random graph when it is a.a.s. Hamiltonian. 
Several recent results about Hamiltonicity ([22], [26], [34], [11]) can be used to show 
fairly easily that G(n,p) with p = p(n) above the threshold for Hamiltonicity contains 
typically many, or even exponentially many Hamilton cycles. Here we aim however for 
(relatively) accurate bounds. 

Using linearity of expectation we immediately see that the expected value of the 
number of Hamilton cycles in G(n,p) is ^ n ^' p n . As the common intuition for random 
graphs may suggest, we expect the random variable to be concentrated around its 
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mean, perhaps after some normalization (it is easy to see that the above expressions 
for the expectation become exponentially large in n already for p inverse linear in n). 

The reality appears to confirm this intuition - to a certain extent. Denoting by X 
the number of Hamilton cycles in G{n,p), we immediately obtain X < (^) a.a.s. 
by Markov's inequality. Janson [20] considered the distribution of X for p = Q (y/n) 
and proved that X is log-normal distributed, implying that X = (^r) n (l + o(l)) n 
a.a.s. It is instructive to observe that assuming p = o(l), the distribution of X is in 
fact concentrated way below its expectation, in particular implying that XfE(X) A 0. 
For random graphs of density p = o (n^ 1 / 2 ) not much appears to be known about the 
asymptotic behavior of the number of Hamilton cycles in corresponding random graphs. 
We nevertheless mention the result of Cooper and Frieze [8], who proved that in the 
random graph process typically at the very moment the minimum degree becomes two, 
not only the graph is Hamiltonian but it has (log n)^ 1 ~ o( - 1 ^ n Hamilton cycles. 

Our main result is the following theorem, which can be interpreted as an extension 
of Janson's results [20] to the full range of p{n). 

Theorem 1. Let G ~ G(n,p) with p > lnn + lnl ^ n + aj ( 1 ) _ Then the number of Hamilton 
cycles is n\p n (l — o(l)) n a.a.s. 

Improving the main result of [8], we also show the following statement. 

Theorem 2. In the random graph process, at the very moment the minimum degree 
becomes two, the number of Hamilton cycles becomes (lnn/e)"(l — o(l)) n a.a.s. 

We continue with a short overview of related results for other models of random 
and pseudorandom graphs. For the model G(n, M) of random graphs with n vertices 
and M edges, notice the result of Janson [20] showing in particular that for the regime 
n 3/2 ^ m <- 0.99 C^), the number of Hamilton cycles is indeed concentrated around 
its expectation. The situation appears to change around M = & (n 3 / 2 ) , where the 
asymptotic distribution becomes log-normal instead. Notice also that the number of 
Hamilton cycles is more concentrated in G(n,M) compared to G(n,p); this is not 
surprising as G(n, M) is obtained from G(n,p) by conditioning on the number of edges 
of G being exactly equal to M, resulting in reducing the variance. 

For the probability space of random regular graphs, it is the opposite case of very 
sparse graphs that is relatively well understood. Janson [19], following the previous 
work of Robinson and Wormald [31], [32], described the asymptotic distribution of the 
number of Hamilton cycles in a random d-regular graph G(n, d) for a constant d > 3. 
The expression obtained is quite complicated, and we will not reproduce it here. For the 
case of growing degree d = d(n), the result of Krivelevich [25] on the number of Hamil- 
ton cycles in (n, d, A)-graphs in addition to known eigenvalue results for G n ^ imply 
an estimation on the number of Hamilton cycles in G n ^ with a superpolylogarithmic 
lower bound on d. 

For an overview of these results as well as of the corresponding results in pseudoran- 
dom settings, we refer the interested reader to [25]. 

1.1. Definitions and notation. The random oriented graph G(n,p) is obtained from 
G(n,p) by randomly giving an orientation to every edge (every of the two possible 
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directions with probability 1/2). Notice that whenever we use the notation G for an 
oriented graph, there exists an underlying non-oriented graph obtained by omitting 
the orientations of the edges of G; it is denoted by G. Making the notation consistent, 
when omitting the vector arrow above an oriented graph, we refer to the underlying 
non-oriented graph. 

Given a graph G, we denote by h{G) the number of Hamilton cycles in G. We call a 
spanning 2-regular subgraph of G a 2-factor. Notice that every connected component 
of a 2-factor is a cycle. We denote by f(G, s) the number of 2-factors in G with exactly 
s cycles. Similarly, a 1-factor of an oriented graph G is a spanning 1-regular subgraph, 
i.e., a spanning subgraph with all in- and outdegrees being exactly one. Analogously, 
the number of 1-factors in G with exactly s cycles is denoted by / ^G, s^j . For the 
purposes of our proofs, we relax the notion of a 2-factor and call a spanning subgraph 
H C G an almost 2-factor of G if H is a collection of vertex-disjoint cycles and at 
most |V(G)|/ ln 2 (|y(G)|) isolated vertices. We denote the number of almost 2-factors 
of G containing exactly s cycles by f'(G, s). Similarly to the notation for non-oriented 
graphs, we call an oriented subgraph H of G an almost 1-factor of G if H is a 1-regular 
oriented graph on at least |V(G)| — |V(G)|/ ln 2 (|F(G)|) vertices. The number of almost 
1-factors of G with exactly s cycles is denoted by /' yG, s^j . 

As usual, in a graph G for a vertex x € V(G) we denote by oIq{x) := \Nq(x)\ its 
degree, i.e., the size of its neighborhood. We denote by 5(G) and respectively A(G) its 
minimum and maximum degrees. For a set S C V(G), we denote by Na(S) the set of 
all vertices outside S having a neighbor in S. Whenever the underlying graph is clear 
from the context we might omit the graph from the index. Similarly, in an oriented 



graph G for a vertex x S V ( G) we call d i q(x 



yeV(G) : yxeE(G 



indegree of x and d Qut ^(x) 



eV G : xzeElG 



the 



the outdegree of x. We 



denote by <5j n yGj , Aj n y3j , 5 ou t \ G) , and A out \ GJ the minimum and maximum in- 
and outdegrees of G. 

In a graph G for two sets A, B C V(G) we denote by ec(A, B) the number of edges 
incident with both sets. In an oriented graph G, for two sets A, B C V(G) the notation 
e^(A, B) stands for the number of edges going from a vertex in A to a vertex in B. We 
write ec(A) := ec(A,A) and e^(^4) := eg(A,A) for short. Similarly to the degrees, 
whenever the underlying graph is clear from the context we might omit the graph from 
the index. 

To simplify the presentation, we omit all floor and ceiling signs whenever these are 
not crucial. Whenever we have a graph on n vertices, we suppose its vertex set to be 

w. 



1.2. Outline of the proofs. In Section 2, the lower bounds for Theorems 1 and 2 are 
proven in the following steps. 

• In Lemma 3 we show using the permanent of the incidence matrix that under 
certain pseudorandom conditions, an oriented graph contains sufficiently many 
oriented 1-factors. 
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• In Lemma 5 we prove that the random oriented graph G(n,p) a.a.s. contains 
a large subgraph with all in- and outdegrees being concentrated around the 
expected value. This subgraph then satisfies one of the conditions of Lemma 3. 

• In Lemma 6 we show that the random graph G(n,p) contains many almost 2- 
factors a.a.s. In the proof, we orient the edges of G(n,p) randomly and apply 
Lemma 3 to the subgraph with almost equal degrees whose existence is guaran- 
teed by Lemma 5 a.a.s. 

• In Lemma 7 we prove that most of these almost 2-factors have few cycles a.a.s. 

• We then call a graph p-expander if it satisfies certain expansion properties and 
show in Lemma 9 that in the random graph process, the graph G(n,p) has these 
properties in a strong way. 

• Lemma 10 shows that in any graph having the p-expander properties and min- 
imum degree 2, for any path Pq and its endpoint v\ many other endpoints can 
be created by a small number of rotations with fixed endpoint v±. 

• Lemma 11 contains the main technical statement of the paper. Its states that 
in a graph satisfying certain pseudorandom conditions, for almost every almost 
2-factor F with few components, there exists a Hamilton cycle with a small 
Hamming distance from F. The proof is a straightforward use of Lemma 10. 

• The proofs of Theorems 1 and 2 are completed with a double counting argument. 
On the one hand, by Lemma 7 there exist many almost 2-factors with few cycles 
a.a.s. Furthermore, for each of these almost 2-factors there exists a Hamilton 
cycle with small Hamming distance from it a.a.s. by Lemma 11. On the other 
hand, for each Hamilton cycle, there are not many almost 2-factors with few 
cycles having a small Hamming distance from it. Hence, the number of Hamilton 
cycles is strongly related to the number of almost 2-factors with few cycles, 
finishing the proof. 



2. The proofs 
Let G ~ G(n,p). Since E(/i(G)) = (n - l)!p n /2, we obtain 

h(G) < In n(n - l)\p n /2 < 

a.a.s., using just Markov's inequality. Thus, for the remainder of the section we are 
only interested in the lower bound on the typical number of Hamilton cycles in the 
random graph. 

We know from [21] and using e.g. the results from [22] and [26] that in G ~ G(n,p) 

there are at least ( [(?(( ; )/2j n! 2 -f actors s. We now want to give an a.a.s. lower 

bound on the number of 2-factors in G, and we want to do it within a multiplicative error 
term of at most 2°^ from the "truth" , basically deleting the 2 from the denominator 
in the above expression in the case p S> lnra/ra, and replacing the term [5(G) /2\ by 
asymptotically np. 

We first prove a pseudo-random technical statement that will give us the desired 
inequality once we show that G (or a large subgraph of it) satisfies the pseudo-random 
conditions. The proof is based on the permanent method as used in [12]. 
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Lemma 3. Let n be sufficiently large and r 3> In Inn, and let G be an oriented graph 
on n vertices satisfying the following (pseudo-random) conditions: 

• 5 in (G),5 ou t(G),A in (G),A out (G) £ (r - 4r/ In In n,r + 4r/lnlnn) 

• for any two sets A,Bc V(G) of size at most \A\, \B\ < 0.6n, there are at most 
0.8ry/\A\\B\ edges going from A to B. 

Then G contains at least ( r 100r / lnlnw N j oriented 1- factors. 



Proof Create an auxiliary bipartite graph G' from G in the following way: take 
two copies X and Y of the vertex set [n] by doubling each vertex v G [n] into v x £ X 
and vy £ Y. We put a (non-oriented) edge uv £ E(G') between vertices ux £ X and 
v y £ y if uv £ E(G) is an edge oriented from ti to u in C We observe a one-to-one 
correspondence between oriented 1-factors in G and perfect matchings in G' . 

In order to use the permanent to obtain a lower bound on the number of perfect 
matchings of G', we need a (large) spanning regular subgraph of G' . Its existence is 
guaranteed by the following claim. 

Claim 4. G' contains a spanning regular subgraph G" with regularity at least d = 
r — lOOr/lnlnn. 

Proof Applying the Ore-Ryser theorem [29] we see that the statement of the 
claim is true provided that for every Y' C Y we have 

d\Y'\ < mm{d,e G/ (x,Y')}. 

Suppose to the contrary that this conidition does not hold, i.e., there exists aY'CY 
s.t. 

d\Y'\ > mm{d,e G ,(x,Y')}. 

We examine the number of edges incident to Y' that can be deleted from G' without 
disturbing the right hand side of the above inequality. Formally, we denote it by 
c = ^2 x& x m ax{0, eQi{x, Y') — d}. Notice that 

c = e(X, Y') - Y min{d, e G > (x, Y')} > e(X, Y') - d\Y'\ 
xex 

as supposed above. 

Since (r -4r/lnlnn)|Y'| < 5{G')\Y'\ < e(X,Y') < d\Y'\ + c, we obtain 

c > rn — y ■ 

In Inn 

On the other hand, denoting by X' the set of vertices that have at least d neighbors in 
Y', and noticing that A(G') < r + 4r/lnlnn, we obtain 

<r 104r iv/i 
m In n 

Hence, 

(1) > 0.9|Y'|. 
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Notice that by the choices of Y' and X' , we have 

(2) d\Y'\ > J2mm{d,e GI (x,Y')} = d\X'\ + e(X \ X' , Y'). 

For the number of edges between Y\Y' and X \ X' we see that 

(r + 4r/ In Inn) | Y \ Y'\ > e{X \ X' , Y \ Y') = e(X \ X',Y) — e(X \ X' , Y') 

(2) 

> 5{G')\X\X'\-d{\Y'\-\X'\) 

(3) > — — \X\X'\ + (r-100r/lnlnn)|Y\Y / |, 

In Inn 

leading to 

(4) \X\X'\ <1.1\Y\Y'\. 
Furthermore, notice that by (2) it holds that 

(5) \X'\ < \Y'\. 

We prove the claim by case analysis. 

• If \Y'\ < n/2, we obtain for the number of edges between X' and Y' 
Choice of x' , (l) 

\ J\ V' I ^ 



e(X',Y') > d\X'\ > 0.9dy/\X'\\Y'\ > 0.8r^/\X'\\Y'\, 

contradicting the second condition of the lemma. 
• If \Y'\ > n/2, then again by the definition of X' we obtain e(X' ,Y') > d\X'\ 
leading to 

(2) (4) (5) 

e(X \ X', Y ) < d(\Y'\ - \X'\) = d(\X \ X'\ - \Y \ Y'|) < 0.1d\Y \ Y'\ < 0.1d\X \ X'\ 
Thus, using the fact that 5(G') > d, we see that 



e(X \X',Y\ V) > 0.9d\X \ X'\ > 0.8ry/\X\X'\ ■ \Y\Y'\, 

again contradicting the same condition of the lemma, since now both X \ X' 
and Y \ Y' have size less than 0.6n by (4). 

□ 

We observe that the number of perfect matchings in G" equals the permanent of the 
incidence matrix of G" . Hence the result of Egorychev [9] and Falikman [10] on the 
conjecture of van der Waerden implies that the number of perfect matchings in G" is 
at least d n n\/n n > (f) n . □ 

In order to use Lemma 3, we first prove the a.a.s. existence of a large subgraph of 
G(n,p) satisfying the degree-conditions of Lemma 3 a.a.s. 

Lemma 5. Let G ~ G(n,p) with p > In njn. Then there exists a set V C [n] 
of at least n — n/ln 2 n vertices of G such that the graph C := G[V] satisfies 

5 m (C),5 out (C),A in (C),A out (C) G ( np^nv/Xnlnn n P+ n P ,lnlnn \ ^ 
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Proof We observe using Chernoff's inequality (see e.g. Corollary A. 1.14 in [2]) 
that for p 3> In n(ln In n) 2 /n the statement holds for V' = [n] a.a.s. Hence, from now 
on we assume np = 0(lnn(lnlnn) 2 ). 

Let L be the set of all vertices whose in- or outdegree is at most np np / lnlnn _|_ \_ p Q r 
every y G [n], we can estimate using Chernoff's inequality (see e.g. Corollary A. 1.14 
in [2]) 

Pr(y G L) = exp (-0 (lnn/(lnlnn) 2 )) . 
Thus, by Markov's inequality we obtain 

\L\ < Inn • E(|L|) = Inn • (n - 1) exp (-Q (lnra/(lnlnn) 2 )) 
(6) = raexp (— (lnn/(lnlnn) 2 )) 

a.a.s. 

Fix an arbitrary vertex x G [n]. We denote 

f r . . . , , . . np — np/lnlnn , . . np—npl In Inn 
L x = |y G [n] \ {x} : d m6 _ x {y) < JL 1 or d out6 _ x (y) < JL 1 

Notice that L x C L, and thus (6) bounds \L X \ as well. 

Since for every y G [n] \ {x} the events "rry G E(G)" and "y G L^" are independent, 
we obtain using Chernoff's inequality again (see e.g. Theorem A. 1.12 in [2]) 



Pr 



(\N G (x) n L x \ > 2 ^ ) I (|£x| = nexp (-0 (lnn/(lnlnn) 2 ))) 



(7) < exp (~ ^f— (lnn/(lnlnn) 2 )J = o(l/n). 



Similarly, we let i? be the set of all vertices whose in- or outdegree is at least 

np—npl In Inn -, l l a • 

— — ^ 1 and obtain 

(8) \R\ < nexp (-0 (lnn/(lnlnn) 2 )) 



a.a.s. 

We define analogously 



J ?ip + np/lnlnn np + np/lnlnrt 

A, = jy £ [n] \ {*} : d in d _ x (y) > 1 or d Qut ^ x {y) > 



and observe analogously to (7) that R x C R and 
(9) 

Pr \(\Na(x)nRr\ > „ 

2 In In n 



|A^g(x) n ^1 > 777^) I (1^1 = ™exp (-0 (lnn/(lnlnn) 2 )))j = o(l/n) 



We denote by V the set of all vertices from [n] whose in- an outdegrees in G lie in 



np—npj In Inn np+np/ In Inn 
2 ' 2 



Notice that [n] \ V C L X L)R X for every x G [n]. Hence, we 
see that | V'| > n — a.a.s. by (6) and (8). Furthermore, from (7) and (9) we obtain 
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that all in- an outdegrees in G[V] lie in (^ np 3n P/ lnlnn ; np+np/ in in n ^ a a s ^ completing 
the proof of the lemma. □ 

From now on, whenever we have n and p chosen, we denote 

d = d(n, p) = np — lOOnp/ In In n . 

In the following lemma, we show that the random graph contains a.a.s. many 2- 
f actors. 

Lemma 6. The random graph G ~ G(n,p) with p > lnn/n satisfies 

^ 2 s f'(G,s) > d~ n ' la2n {d/e) n 
s€[n/3] 

a.a.s. 

Proof In order to use Lemma 3, we orient G at random to obtain G (as always 
in this paper, for every edge each of the two possible orientations gets probability 1/2 
independently of the choices of other edges). 

First we show that the second condition of Lemma 3 holds a.a.s. for G with the 
intuitive choice r = np/2. Since the maximum degree of G is at most 3np a.a.s. (see 
e.g. [5]), we obtain that in G a.a.s. for any two sets A and B with \A\ > 100|-B| the 
number of edges between them is at most 3np\B\ < 0. Anp^/\A\ \B\ . Hence, e ( =i(A, B) < 

0Anpy/\A\\B\ and eg(B, A) < 0Anpy/\A\\B\. Thus, we are left with the case of sets 
A and B of sizes \ A\ < W0\B\ and \B\ < W0\A\. 

For small disjoint sets, we obtain using Chernoff's bound (see e.g. Theorem A. 1.12 
in [2]) 



/ n 2 

Pr 3A', B' C [n], A' n B' = 0, \A'\\B'\ < — , \A\ < 100151 < 10 4 |A| : eMA',B') > OAnpJlA 7 ] 

\ In In n u 



a,b=o(n), a=0(b) 



n \ / n — a 



exp —0 npyab In 



npV ah 
pab 



< 



E 



t)°(t) ^(-"(""'""'W; 



a,b=o(n) 

< Y ex p ( aln (— ) + &ln (-^) _ ° ( aln (-) lnn ) _ ^ ( &ln (^-) lnn 

a,b=o(n) 

= o(l). 



nfftlnnlnin 



n 
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Similarly, for large disjoint sets we obtain using Chernoff's bound (see e.g. Corol- 
lary A. 1.14 in [2]) 

/ n 2 

Pr 3A',B' C fnl, A' n B' = 0, \A'\\B'\ > — , \A\,\B\ < 0.6?i : eMA'.B 1 ) > 0AnpJ\A'\\B' 

V In In n u 



s E (;)("; a )-P(-n(*)) 



2 

a,b<n, ab>r^ — 

' — ' In In n 

( ^ ( nlnn \ \ 

24 exp (" f! (^)) =0(1) ' 

Hence, a.a.s. for every pair of disjoint sets A' and B' , the number of edges going from 
A' to B' satisfies 

(10) e & {A', B') < 0Anpy/\A'\\B'\. 
Analogously, we see that a.a.s. for every M C [n] of size at most 0.6n, 

(11) e G (M) < 0Anp\M\. 

Thus, a.a.s. for every A, B C [n] of size \A\, \B\ < 0.6n, the number of edges going 
from A to B in G is bounded by 

(io), (ii) . 

eg(A,B) =e a (A\B,B\A) + e G (AnB) < 0Anpy/\A \ B\ \B \ A\ + 0Anp\ A n B\ 



< 0Anpy/\A\\B\, 

establishing that the second condition of Lemma 3 holds a.a.s. for every subgraph of 
G. 

Hence, by Lemma 5 the graph G ~ G(n,p) a.a.s. is such that for a random orienta- 
tion G, there a.a.s. exists a vertex set V' C [n] of size at least n — nj In 2 n such that the 
induced subgraph G[V] satisfies the conditions of Lemma 3 with r = np/2. Applying 
Lemma 3 to this induced subgraph, we obtain 



se[n/3] 



a.a.s. Thus, we obtain 

n—n/ In 2 n 



sS[n/3] 

a.a.s., where the expectation is taken over the random choice of orienting the edges of 
G, the process creating G from G. 

On the other hand, when we orient the edges, an almost 2-factor of G with exactly 
s cycles becomes an almost 1-factor of G with probability at most 2 in 2 n ~ n+s , implying 

£ 2Y(G, S )> £ 2 re -I^E (/'((?,. 

se[n/3] se[n/3] 



10 



R. GLEBOV AND M. KRIVELEVICH 



Putting these two facts together, we obtain 

(r\ n—n/ In 2 n 

a.a.s., completing the proof of the lemma. □ 

We show now that there are typically many almost 2-factors in G with a small 
number of cycles. We denote 

n 



s* = s*(n) 



In nVhi In n 

Lemma 7. For every p > \nn/n, the random graph G ~ G(n,p) satisfies 

s* 

Y,f'{G,s)>(np/e) n {l-o(l)Y 

s=l 

a.a.s. 

Proof By Lemma 6 we know that E se [n/3] 2S f'( G , s ) > d~ n / ln2 n (d/e) n a.a.s. 

We show now that the contribution of almost 2-factors with too many cycles is 
negligible. We use the estimate (5) of [21]: in the random graph H ~ G(n',p), for 
every s > Inn', 

E(/ ( H..)) <'"'- : >! (in "' r ' p ~. 

Ul ' " - (s- 1)!2 S 

We obtain 

n/3 , . n/3 , . s £ 

n \ ^ n\ (InnJ p 



X;E(2-m-))< E HE 



l<n/\n 2 n «=s* 



. n/\n 2 n n / 3 

n \ 1 v-^ ( s 



£ 

9=S* 

(d/e) n e° (n/lDlnn) ( — 



v J ^— ' Velnn 

/ S=S* 

* \ —s 



e In n 



o(d- n / ln2n (d/e) n ) 



Hence, using this estimate together with Markov's inequality, we see that the number 
of almost 2-factors of G with at most s* cycles is 

s* 

jh f'(G, s) > l -2~ s 'd- n '^ n {d/e) n = (np/e) n (1 - o(l))» 
a.a.s. □ 



The next technicality we need to prove in order to be ready to prove the main 
theorem is the expansion of G(n,p). 
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To collect all but one expansion properties that we need, we make the following 
definition. 

Definition 8. We call a graph G with the vertex set [n] a p-expander, if there exists a 
set D C [n] such that G and D satisfy the following properties: 

• \D\ < n om . 

• The graph G does not contain a non-empty path of length at most 3 ^"" n such 
that both of its (possibly identical) endpoints lie in D. 

• For every set S C [n] \ D of size \S\ < -, its external neighborhood satisfies 
\N(S)\>^\S\. 

The following lemma shows that these properties are pseudo-random. 

Lemma 9. Consider the two-round expansion of the random graph and fix G ~ G(n,p) 
with lnn/n < p < 1 — 21nlnra/ra and GCG~ G(n,p) with p = p + 21nlnn/n. Then 
it is a.a.s. true that every graph G' satisfying G C G' C G is a p-expander. 

Proof We first expose G and fix D = {v 6 [n] : d G {v) < np/100} to be the 
set of all vertices of G with degree less than rap/100 in G. Since for a fixed set D the 
second property is decreasing and the third property is increasing, it suffices to prove 
the second statement for G and the third statement for G. 

The first property is satisfied by Claim 4.3 from [4] a.a.s. The second property can 
be proven to hold in G a.a.s. similarly to Claim 4.4 from [4] (there it is proven to hold 
for G a.a.s.) 

For the third property, assume to the contrary that there exists a set S C [ra]\D of size 
at most |5| < - such that its external neighborhood in G satisfies \Ng(S)\ < ^q\S\. 
By the definition of D, the number of edges incident to S in G is 

e G (S,N G (S)US) > \S\np/200. 

But Chernoff's inequality (see e.g Theorem A. 1.12 in [2]) tells us that 

Pr (3A, B<Z[n], \A\ < \ \B\ < ^\A\ : e G (A, B U A) > \A\np/200 



v , n \fn-\A\\( B(\e G (A,BuA)\) \ 
V|A|yV 1^1 )\ |^|rap/200 J 



< 

A,Bc[n],\A\<l/p, \B\<\A\np/1000 



\A\np/200 



A,Sc[n],|A|<l/p,|B|<|A|np/1000 Vl 17 Vl 17 V 



ra \ / n \ f 200e\A\\AU B\p 



\A\np 



\A\np /200 



( 3ap \ 

j I anp_ lie) 



anp/200 



< ^ a> n P\ n J \ anp 
a<l/p 

< an P ( -T- ) = 



a<l/p 



^3apy np/40 ° 
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providing that the third property holds in G a.a.s. □ 

The proof of the next lemma is based on the ingenious rotation-extension technique, 
developed by Posa [30], and applied later in a multitude of papers on Hamiltonicity, 
mostly of random or pseudorandom graphs (see for example [7], [12], [23], [27]). 

Let G be a graph and let p> = (v±, i>2, • • • , v q ) be a path in G. If 1 < i < q — 2 and 
(v q , Vi) is an edge of G, then there exists a path P' = (v\V2 ■ ■ ■ ViV q v q -\ . . . i>j+i) in G 
with the same set of vertices. The path P' is called a rotation of Po with fixed endpoint 
v\ and pivot Uj. The edge is called the broken edge of the rotation. We say 

that the segment t> j+i . . . v q of Po is reversed in P 1 . In case the new endpoint t> j+i has 
a neighbor Vj such that j £ {i,i + 2}, then we can rotate P' further to obtain more 
paths of the same length. We will use rotations together with the expansion properties 
from Lemma 9 and the necessary minimum degree condition to find a path on the same 
vertex set as Po with large rotation endpoint sets. 

The next lemma shows that in any graph having the p-expander property and min- 
imum degree 2, for any path Po and its endpoint v±, after a small number of rotations 
with fixed endpoint v\, we either create many other endpoints or extend the path. Its 
proof has certain similarities to the proofs of Lemma 8 from [15] and of Claim 2.2 
from [17]. 

Lemma 10. Let n be a sufficiently large integer and G be an n-vertex p- expander with 
minimum degree 5(G) > 2 and np > Inn. Let Pq be a viw-path in G. Denote by 
B{v\) C V(Pq) the set of all vertices v € V for which there is a v\v-path on the vertex 
set V(Pq) which can be obtained from Pq by at most 3 ^7^3 rotations with fixed endpoint 
V\. Then B(v\) satisfies one of the following properties: 

• there exists a vertex v 6 B(v\) with a neighbor outside V(Pq), or 

• \B(vi)\ > n/3000. 

Proof Assume that B(v\) does not have the first property (i.e., for every v £ 
B(vt) it holds that N(v) C V(P )). 

Let to be the smallest integer such that (3500)*° 1 > |; note that to < 2 lr |° n ^ • Since 
G is a p-expander, there is a corresponding vertex set D as in Definition 8. 

At the first step, we find a neighbor u D U N(D) of w that is not a neighbor of w 
along Po. Its existence is guaranteed since w has at least two neighbors along Po, and 
by the second p-expansion property, at most one of them can have a neighbor in D. 
We rotate the initial path Po with pivot u and call the resulting path P' = (v\, . . . , v q ). 
Notice that this way, v q is guaranteed not to belong to D. 

We construct a sequence of sets So, . . . ,St Q B(v{) \ D C V(Po) \ {^1} of vertices, 
such that for every < t < to and every v € St, v is the endpoint of a path which can 
be obtained from P' by a sequence of t rotations with fixed endpoint v\ , such that for 
every < i < t, the non-f i-endpoint of the path after the ith rotation is contained in 
Si. Moreover, \S t \ = (3^5)* for every t < t - 2, \S k) -i\ = |, and |5 to | > 3^0. 

We construct these sets by induction on t. For t = 0, one can choose So = {v q } and 
all requirements are trivially satisfied. 
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Let now t be an integer with < t < to — 1 and assume that the sets So, ... , St—i 
with the appropriate properties have already been constructed. We will now construct 
St- Let 

t-i 

Vi G N(St-i) : Vi, v i+ i g" [J Sj U D 

j=0 

be the set of potential pivots for the tth rotation, and notice that T C V{Po) due to our 
assumption, since T C N(St-i) and St—i C B{v\). Assume now that vi £ T , y £ St-i 
and (vi,y) £ E{G). Then, by the induction hypothesis, a uiy-path Q can be obtained 
from P' by t — 1 rotations such that after the jth rotation, the non-vi-endpoint is in 
Sj for every < j < t — 1. Each such rotation breaks an edge which is incident with 
the new endpoint, obtained in that rotation. Since Vi-\,Vi,Vi + \ are not endpoints after 
any of these t — 1 rotations, both edges (vi-\,Vi) and (vi, fi+i) of the original path P' 
must be unbroken and thus must be present in Q. 

Hence, rotating Q with pivot Vi will make either Vi-\ or an endpoint (which of 
the two, depends on whether the unbroken segment Vi-iViVi+i is reversed or not after 
the first t — 1 rotations). Assume without loss of generality that the endpoint is 
We add V\-\ to the set St of new endpoints and say that Vi placed Vi-\ in St- The only 
other vertex that can place in St is Uj_2 (if it exists). 

Observe now that if t < 0.1 Inn/ In Inn, the distance between any vertex from 
and v q is at most 2t — 2 < 0.2 Inn/ In Inn by the way the sets were constructed. Hence, 
between any two vertices from A r (5t_i)U N(N(St-i)), there is a path of length at most 
0.51nn/lnlnn. Thus at most one vertex from D can be in N(St-i) U N(N(St—i))- 
On the other hand, it t > 0.1 In nj In Inn, then \D\ < n om = o(\S t -i\) = o(\N (S t -i)\) - 
Thus, in both cases \D n (N(S t -i) U N(N(S t -i)))\ = o(\N(S t -i)\). 

Combining all this information together, we obtain 

I & I > -\T\ 
2 1 

> \ (|JV(5t-0l - 3(1 + | Si | + • • • + + \D n (iV(5t_i) U #(#(&_!)))!)) 
^ np \* 



3000 



Clearly we can delete arbitrary elements of St to obtain St of size 

(3S00) if * < *o-2 

and of size - if t = to — 1. So the proof of the induction step is complete and we have 
constructed the sets So, • • • > St -i- 

To construct St we use the same technique as above, only the calculations are 
slightly different. 

\S t0 \>l\T\ 

> 1 (|iV(S t0 _i)| - 3(1 + | Si | + . . . + |S 4o _ 2 | + |S t0 _i| + \D n (N(S t -i) U N(N(St-i)))\)) 

> n/3000. 



11 



R. GLEBOV AND M. KRIVELEVICH 



The set St ■= St is by construction a subset of B(vi), and the number of rotations 
needed to make any of its vertices an endpoint of the current path is at most to + 1, 
concluding the proof of the lemma. □ 

The proof of the following lemma relies on the final part of the proof of Theorem 1 
from [25] and uses Lemma 10. It shows that under certain pseudorandom conditions 
in a graph G for every almost 2-factor, after adding few random edges, there exists a 
Hamilton cycle within a small Hamming distance from it a.a.s. 

Lemma 11. Let G be a connected n-vertex p-expander with minimum degree 2 and S 
be a set of vertices of G of size \S\ = o(n) such that there exist at least n non-edges in 
G not incident to S. Let F be an almost 2-factor of G with at most s* cycles. Choose 
n non-edges e±,...,e n of G i.a.r. under the condition that none of them is incident to 
S and denote by G' the (random) graph obtained from G by turning them into edges. 
Then, if it is a.a.s. true that every graph G C G C G' is a p-expander, then G' a.a.s. 
contains a Hamilton cycle H with Hamming distance at most 17s* lnn/ln(np) from F. 

Proof Fix an arbitrary component CCf, Since G is connected, there exists 
an edge in G connecting a vertex v G V(C) and y V(C) - unless of course C is 
already Hamiltonian. We denote by C" the component of y in F. Opening C up by 
deleting an edge of C incident to v (no need to do so if C is just one isolated vertex), 
we get a path P. We append the edge vy to P, go through it to C', and if C' is a 
cycle, then we open it up by deleting an edge of C' incident to y to get a longer path P' 
and repeat the argument. If at some point there are no edges between the endpoints 
of the current path P" and other components from F, then we can fix one endpoint x 
of P" and rotate P" using Lemma 10 to extend it outside or to obtain a set B{x) of 
size at least |J5(cc)] > n/3000 of potential other endpoints. For every vertex z G B(x), 
we can rotate the resulting path fixing z as one endpoint to obtain a set A(z) of size 
at least |-A(z)| > n/3000 of potential other endpoints or to extend the path outside. If 
the path still cannot be extended outside and we can still not close it to a cycle, we 
have a set E' of at least 10 _8 n 2 non-edges of G not incident to S, so that turning any 
of them into an edge would close the path to a cycle. We add pairs e±, e^i ■ ■ ■ to E(G), 
until one of them falls inside E' . Notice that for every i € [n], the pair ei falls into E' 
with at least some constant positive probability. This means that considering events 
"ej S E'" , every event has probability 6(1) regardless of the previous events. Notice 
that in a successful round, the number of components gets reduced or a Hamilton 
cycle is created, since the edge that appeared in E' closed the path into a cycle or 
extended the path directly. To reduce the number of components by one, we do at 
most lnn/ln(np) rotations by Lemma 10, therefore increasing the Hamming distance 
from F by at most 4 + 12 In n/ ln(np). Since it is enough to have s* +n/ In 2 n successful 
events to obtain a Hamilton cycle, the expected number of needed turns of non-edges 
into edges is at most 0(1) • (s* + n/ln 2 n) = o(n). Hence the n additional edges suffice 
to create a Hamilton cycle H from F by Markov's inequality a.a.s., replacing at most 
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81nn/ln(np) edges for every component of F. Thus, the Hamming distance between 

; In n , 
ln(np) 



F and H is at most 2 • 8 r ^ T (s* + nj In 2 n) < 17s* Inn/ In (np). □ 



We are now ready to prove Theorem 1. 

Proof Notice that only the lower bound is of interest for us. We expose G in 
two rounds. 

We choose a function p\ = p\{n) such that lnn + lnl " n + aJ ( 1 ) < p 1 < p — ^p-, p\ < 
l?21nlnn/n, and p± = (1 — o(l))p. In the first round, we expose G\ ~ G(n,pi). We 
determine D := {v E [n] : dd(v) < npi/100}. 

In the second round, we expose the binomial random graph G2 by including every 
edge from K n \ G\ into E(G2) with probability P2 ■= fEfj- Since (1 — Pi)(p2) = 1 — P, 
we obtain a graph G := G\ U G2 ~ G(n,p). We know by Lemma 9 that a.a.s. every 
graph between G\ and G including them both is a pi-expander. Furthermore, notice 
that the expected number of edges in G2 is ( 7 2 l )p2 > GjKl — Pi) = w ( n )> hence a.a.s. 
at least n additional random edges appeared in the second round of expansion by 
Markov's inequality. Since these edges were chosen i.a.r., and G\ is a.a.s. connected 
with minimum degree at least 2 (see e.g. [5]), the conditions of Lemma 11 are satisfied 
for G\ and the first n edges exposed in the second round with 5 = 0. 

Now, we put all we know together: 

• By Lemma 7 we obtain X^s=i f'i^i, s) > (npi/e) n (1 — o(l)) n a.a.s. 

• For every almost 2-factor F of G\ with at most s* cycles, there a.a.s. exists a 
Hamilton cycle in G with Hamming distance at most k := 17s* lnn/ln(npi) = 



Yin 



o(n) from F by Lemma 11. 



ln(npi ) %/lnlnn 

• On the other hand, for every Hamilton cycle H in G, to obtain an almost 2- 
factor of G\ of distance at most k from H, we can first delete at most k edges of 
H, thus obtaining a collection of at most k paths. These paths should then be 
tailored into an almost 2-factor, and the choices here are for each of the at most 
2k endpoints of the paths to be connected to one of its A(Gi) neighbors in G\ 
or to stay isolated. Thus, there are at most (£)(A(Gi) + l) 2fe almost 2-factors 
of G with Hamming distance at most k from H. 

• Hence, by double counting almost 2-factors of G with at most s* cycles, we 
obtain 

h(G) >- feST& £ (n! " /e) " A^Z))» = (npleT (1 " 0<1)) " 

a.a.s. □ 



To strengthen the result of Cooper and Frieze [8], we now prove Theorem 2. 

Proof The proof goes along the argument of Theorem 1, but now we expose 
the graph in three rounds. We first expose G\ ~ G(n, lnn/n) and fix the set D of 
vertices of degree at most Inn/ 100. Notice that similarly to the argument in the proof 
of Lemma 9, Claim 4.3 from [4] implies that \D\ < n 009 . In the second round of 
exposure, in addition to G\ we expose those edges that are incident to D one by one, 
until the minimum degree becomes two; the resulting graph is called G' . In the third 
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round of exposure we consider the binomial random graph Gi by including every edge 
of K n \ G\ not incident to D with probability p2 := i" n . 

Let us denote by G the graph obtained by stopping the random graph process at 
the moment the minimum degree becomes two. Notice first that since in the random 
graph process 5 (G (n, lnn+0 n 51nlnn )) = 1 a.a.s., we obtain G'uG 2 CG a.a.s., where 
by the union of two graphs with vertex sets [n] we denote the graph on the same 
vertex set where the union is taken over the edge sets. Furthermore, observe that since 
in the random graph process 5 (G (n, lnrt+2 n Inlnn )) > 2 a.a.s., we obtain G'UG 2 C 
G (n, lnn+21nlnn ) a.a.s. (The two statements above can be found e.g. in [5].) Finally, 
G' is connected a.a.s. because of the expansion properties and the fact that the edge 
set between two linearly large sets is not empty (see e.g. [17]). 

Since p2 = ui(l/n) and \D\ = o(n) a.a.s., we obtain \E(G2)\ = uj(n) a.a.s. Further- 
more, these edges are random under the only conditions of being non-edges of 67 1 and 
being not incident to D. Hence, the conditions of Lemma 11 are satisfied for G' and 
the first n edges exposed in the third round. 

Thus, following the lines of the proof of Theorem 1, we obtain the desired estimate. 

□ 



3. Concluding remarks 

In this paper we have proven that for any value of the edge probability p = p(n), for 
which the random graph G G(n,p) is a.a.s. Hamiltonian, the number of Hamilton cycles 
in G is n!p n (l + o(l)) n a.a.s., thus being asymptotically equal to the expected value - up 
to smaller order exponential terms. Of course, it would be very nice to extend Janson's 
result [20] to smaller values of p and to understand more accurately the distribution 
of the number of Hamilton cycles in relatively sparse random graphs. However, given 
that the machinery used in [20] is rather involved, and the result (limiting distribution) 
is somewhat surprising, this will not necessarily be an easy task. 

Our bound on the number of Hamilton cycles in G(n,p) can be used to bound 
the number of perfect matching similarly to [25]. Let m(G) denote the number of 
perfect matchings in the graph G. Since every Hamilton cycle is a union of two perfect 
matchings, we obtain h{G) < ( m( 2 G) ). Hence, for G ~ G(n,p) the a.a.s. lower bound 
on h(G) from Theorem 1 provides the a.a.s. lower bound m(G) > (np/e) n / 2 (l — o(l)) n . 
Since the upper bound is easily obtained from the expected value by Markov's inequality 
similarly to the first paragraph of Section 2, we have m(G) = (np/e) n ^ 2 (l — o(l)) n . The 
corresponding hitting time statement is obtained by a straightforward modification of 
the proof of Theorem 2: In the random graph process, the edge that makes the graph 
connected a.a.s. creates (lnn/e) n / 2 (l — o(l)) n perfect matchings. 
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